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We report the discovery of interstellar mercapto radicals along the sight-line to the submillimeter continuum source W49N. We have 
used the GREAT instrument on SOFIA to observe the 1383 GHz ^11^/2 7 = 5/2 <— 3/2 lambda doublet in the upper sideband of the 
LI receiver. The resultant spectrum reveals SH absorption in material local to W49N, as well as in foreground gas, unassociated with 
W49N, that is located along the sight-line. For the foreground material at velocities in the range 37 - 44 km/s with respect to the 
local standard of rest, we infer a total SH column density ~ 4.6 x 10'^ cm"-, corresponding to an abundance of ~ 7 x 10"' relative 
to H2, and yielding an SH/H2S abundance ratio ~ 0.13. The observed SH/H2S abundance ratio is much smaller than that predicted 
by standard models for the production of SH and H2S in turbulent dissipation regions and shocks, and suggests that the endothermic 
neutral-neutral reaction SH + H2 — » H2S + H must be enhanced along with the ion-neutral reactions believed to produce CH+ and 
SH* in diffuse molecular clouds. 
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1. Introduction 

In the seven decades following the discovery of interstellar CH 
(Swings & Rosenfeld 1937) - the first identified interstellar 
molecule - five or six additional neutral diatomic hydrides have 
been discovered in the interstellar gas: OH (Weinreb 1963), H2 
(Carruthers 1970), HCl (Blake et al. 1985), NH (Meyer & Roth 
1991), HF (Neufeld et al. 1995), and SiH (Schilke et al. 2001; 
a tentative detection). These discoveries were obtained from ob- 
servations over a remarkably wide range of wavelengths, from 
the far-ultraviolet (H2, discovered at 1013 - 1110 A), through 
the near-UV (NH at 3358 A), the optical (CH at 4300 A), the far- 
infrared (HF at 122 yum), the submillimeter (SiH at 478 yum, and 
HCl at 479 i-im), and the radio (OH at 18 cm) spectral regions. In 
addition to these neutral molecules, four hydride molecular ions 
have now been discovered: CH* (Douglas & Herzberg 1941), 
OH+ (Gerin et al. 2010a; Wyrowski et al. 2010), SH+ (Benz et 
al. 2010), and HCF (DeLuca et al. 2012; Gupta et al. 2012), the 
latter three in just the past two years. 

The diatomic hydrides represent the simplest of interstellar 
molecules, and - carefully interpreted in the context of astro- 
chemical models - may provide key information about the in- 
terstellar environment. For example, the OH* abundance yields 
constraints upon the rate of cosmic ray ionization, which initiates 
the chemical network leading to OH* (e.g. Neufeld et al. 2010). 
The CH*, SH* and OH abundances, by contrast, probe the in- 
fluence of shocks and turbulent dissipation, which heat the in- 
terstellar gas and/or lead to ambipolar diffusion; these processes 
may thereby drive endothermic chemical reactions that lead to 



enhanced abundances of CH*, SH* and OH (e.g. Flower, Pineau 
des Forets, & Hartquist 1985; Godard, Falgarone & Pineau des 
Forets 2009) 

The mercapto radical, SH - which, like CH* and SH*, is 
expected to show a strong abundance enhancement in warm re- 
gions that are heated by shocks or turbulent dissipation - has 
been conspicuously absent from the list of previously-detected 
molecules. In its ground rotational state, SH can absorb ra- 
diation in the ^113/2 7 - 5/2 <— 3/2 lambda doublet near 
1383 GHz, a frequency that is completely inaccessible from 
the ground, due to atmospheric absorption, and which - by bad 
fortune - falls right in the gap between Bands 5 and 6 of the 
Herschel Space Observatory'' s Heterodyne Instrument for the 
Far-Infrared (HIFI). Fortunately, this spectral region is now ac- 
cessible to the GREAT (German Receiver for Astronomy at 
Terahertz Frequencies) instrument on SOFIA (the Stratospheric 
Observatory for Infrared Astronomy), and in this Letter we re- 
port the first detection of interstellar SH with the use of that in- 
strument. The observations and data reduction are described in 
§2, and the results presented in §3. In §4 we discuss the im- 
plications of the measured SH column density in the context of 
astrochemical models. 



' GREAT is a development by the MPI fiir Radioastronomie and 
the KOSMA / Universitat zu Koln, in cooperation with the MPI fiir 
Sonnensystemforschung and the DLR Institut fiir Planetenforschung. 
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Table 1. ^Hs/a 7 = 5/2 ^ 3/2 SH lines observed with GREAT 



Transition Frequency A,j Av 

(GHz) (s-') (kms-') 



F 


= 2+ 


^ 1- 


1382.9056 


0.47 X 10- 


-3 


+0.98 


F 


= 3+ 


^ T 


1382.9101 


4.72 X 10" 


-3 




F 


= 2+ 


«- r 


1382.9168 


4.24 X 10- 


-3 


-1.45 


F 


= 2- 


«-2+ 


1383.2365 


0.47 X 10- 


-3 


+ 1.02 


F 


= 3- 


«-2+ 


1383.2412 


4.72 X 10- 


-3 




F 


= T 


^ r 


1383.2478 


4.24 X 10- 


-3 


-1.43 



Velocity shift relative to the strongest hyperfine component of the 
respective doublet 



2. Observations and data reduction 

We observed the ^113/2 7 - 5/2 <— 3/2 transition of SH, a 
lambda doublet for which the strongest hyperfine components 
(F = 3 <- 2) lie at frequencies of 1382.910 and 1383.241 
GHz, in the upper sideband of the LI receiver. A broad atmo- 
spheric ozone feature is present in the lower sideband, at a fre- 
quency of 1379.50 GHz. The LO settings were selected to sep- 
arate this feature as far as possible (in IF frequency) from the 
target SH transition. These observations, with a combined on- 
source integration time of 4.0 min, were carried out on 2011 
September 28th as part of the SOFIA Basic Science Program. 
The telescope beam, of diameter ~ 21" HPBW, was centered on 
W49N at Of = 19h 10m 13.2s, 5 = +09" 06' 12.0"(J2000). The 
observations were performed in dual beam switch mode, with a 
chopper frequency of 1 Hz and the reference positions located 
75" on either side of the source along an East- West axis. The 
AFFTS backend provided 8192 spectral channels at a spacing 
of 183.1 kHz. Thanks to laboratory spectroscopy performed by 
Morino & Kawaguchi (1995) and by Klisch et al. (1996), the SH 
rest frequencies are known to estimated accuracy of < 2 MHz, 
The separation of the lambda doublet corresponds to a velocity 
shift of 71.8 km s"', and each doublet member is further split 
into three hyperfine components; the rest frequencies and spon- 
taneous radiative rates are listed in Table 1, for an assumed SH 
dipole moment of 0.758 D (Meerts & Dynamus 1975). 

The raw data were calibrated to the ("forward beam 
brightness temperature") scale, fitting independently the dry and 
the wet content of the atmospheric emission. Here, the assumed 
forward efliciency was 0.95 and the assumed beam efliciency 
for the LI band was 0.54. The uncertainty in the flux calibration 
is estimated to be ~ 20% (Heyminck et al. 2012). Additional 
data reduction - performed using CLASSEI - entailed second- 
order baseline removal, averaging the data (with a weighting in- 
versely proportional to the square of the r.m.s. noise), and spec- 
tral smoothing to a 3.7 MHz channel spacing. 

In addition to the SOFIA observations of SH that are the 
primary subject of this Letter, we have carried out ancillary ob- 
servations of the 1 10 - loi 168.763 GHz transition of H2S, using 
the IRAM 30 m telescope located in Pico Veleta near Granada 
(Spain) in December 2006 under good weather conditions. We 
used the CI 50 receiver, tuned in single side band, coupled to 
two spectrometers: the high resolution correlator VESPA with 
40 kHz spectral resolution, and a broad band filter bank with 
1 MHz spectral resolution. The data were acquired using the 
wobbler with a frequency of 0.5 Hz, for a total time of 32 min- 
utes, and analyzed with the CLASS software. 



Continuum and Line Analysis Single-dish Software - 
p^//www.iram.fr/IRAMFR/GILDAS, 
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8.0 - 




Frequency (GHz) 

Fig.l. Spectrum of SH ^X\-ip_J = 5/2 ^ 3/2 obtained by 
GREAT toward W49N. Note that because GREAT employs dou- 
ble sideband receivers, the complete absorption of radiation at a 
single frequency will reduce the measured antenna temperature 
to one-half the apparent continuum level. The lambda doubling 
and hyperfine splittings are indicated by the red bars for a com- 
ponent at an LSR velocity of 40 km s"' . 



3. Results 

Figure 1 shows the observed spectrum of SH ^113/2 7 = 5/2 <— 
3/2, with the frequency scale coiTected to the Local Standard 
of Rest (LSR). The double sideband continuum antenna tem- 
perature is rj(cont) = 7.2 K, and the r.m.s noise is 0.12 K in 
a 3.7 MHz channel. Because GREAT employs double sideband 
receivers, the complete absorption of radiation at a single fre- 
quency will reduce the measured antenna temperature to one- 
half the apparent continuum level. 

In Figure 2, the fractional transmission is shown separately 
for each of the lambda doublets (top two panels), with the fre- 
quency scale expressed as Doppler velocities relative to the 
Local Standard of Rest (LSR) for the strongest hyperfine com- 
ponent of each doublet member. The transmission is given by 
2r*/r^(cont)-l, given the assumption that the sideband gain ra- 
tio is unity. Analogous spectra are shown for the 1 10- loi (ground 
state) transition of ortho-H2S, and for the ^113/2/ = 3/2" <— 
^ni/2 7=1/2+ transition of CH (Gerin et al. 2010b), beUeved to 
be a good tracer of H2 (e.g. Sheffer et al. 2008). 

Absorption by SH is clearly detected in the range vlsr ~ 
5 - 20 km s"\ near the systemic velocity of the source. In addi- 
tion, a narrow absorption feature is detected unequivocally near 
i^'LSR ~ 39kms"', a component that is clearly present in the 
spectra of CH, H2S, and many other molecules (e.g. Godard 
et al. 2010; Sonnentrucker et al. 2010). This component arises 
in a foreground cloud unassociated with W49N, which has an 
(kinematically-) estimated Galactocentric distance of ~ 6.7 kpc 
(Godard et al. 2012). 

In contrast to the case of CH, there is an absence of strong SH 
or H2S absorption in the 60 - 65 km s ' range. This behavior is 
similar to that observed for CS (Miyawaki, Hasegawa & Hayashi 
1988) and for the nitrogen hydrides NH, NH2, and NH3 (Persson 
et al. 2012); it may suggest that those molecules - such as SH - 
for which the vlsr ~ 39 km s ' absorption is much stronger than 
that in the 60 - 65kms ' range all originate in material with a 
larger molecular fraction. 

Because of its large spontaneous radiative decay rate (4.7 x 
10"^ s"'), the SH transition we have observed possesses a high 
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Table 2. Column densities for vlsr in the 37 - 44 km s ' 
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Fig. 2 

5/2 < 

3/2 <— 7 = 1/2^). The CH spectrum has been hyperfine 

deconvolved. 



Ratio of flux to continuum flux, for SH (^03/2 J - 
■ 3/2), ortho-H2S (lu, ^ loi), and CH CU3/2J = 



critical densit}0 ( i 10^ cm"^, for an assumed collisional deexci- 
tation rate ~ 10"'*'cm''s"', similar to that computed by Off'er et 
al. 1994 for the analogous transition of OH). Thus, in the fore- 
ground material unassociated with W49N, we expect that SH 
will be almost entirely in its ground rotational state, ^113/2/ - 
3/2. In that case, the absorption optical depth, integrated over 
velocity and summed over the the six components listed in Table 
1, is given by 



TC/V 



X 10"''*Ar(SH)cm2kms"'. 



(1) 



(For the 5-20 km s"' velocity interval, equation (1) would likely 
yield an underestimate of the SH column density, because - for 
gas associated with the dense W49N cloud itself - there is likely 
to be a significant SH population in excited rotational states.) 

In the present studjO we confine our attention to the nar- 
row absorption feature appearing near vlsr - 39kms"'. In 
Table 2, we present estimates of the SH column density in the 



Here, "critical density" is defined as the gas density at which the 
coUisional deexcitation rate equals the spontaneous radiative decay rate. 

A future study, requiring detailed modeling of the excitation of SH 
in W49N, will be needed to obtain an estimate of the SH+ column den- 
sity in the source itself. 



Molecule 


Column 


Abundance'' 






density 


relative to H2 






(cm-2) 






SH 


4.6 X 10" 


6.9 X 10-" 


Present work 


ortho-H2S 


2.6 X 10'^ 


3.9 X lO-'* 


Present work 


SH+ 


2.6 X lO'^ 


3.9 X 10-" 


Godard et al. 2012 


para-H20 


1.5 X 10" 


2.3 X 10-** 


Sonnentrucker et al. 2010 


CH 


5.8 X 10" 


9.0 X lO-'* 


Gerin et al. 2012 



For an assumed H2 column of 6.6 x 10^" cm ^ (Gerin et al. 2012) 

39kms-' absorbing cloud, along with analogous results ob- 
tained for ortho-H2S from IRAM 30m observations and for SH^, 
para-H20 and CH from Herschel/HIFI observations. Given an 
estimated H2 column density of 6.6 x lO^^cm"^ for this cloud 
(Godard et al. 2012), the observed SH column density implies an 
SH/H2 abundance ratio ~ 7 x 10"^, corresponding to ~ 0.003% 
of the solar abundance of elemental sulfur. The SH/SH^ ratio 
is ~ 1.8. For an assumed H2S ortho-to-para of 3 (the value ex- 
pected in equilibrium), the SH/H2S ratio is 0.13, a value that is 
signficantly smaller than the OH/H2O ratio -1.0 observed by 
SOFIA for this absorbing cloud (Wiesemeyer et al. 2012) 

4. Discussion 

Sulfur is unusual among the abundant elements in having a set of 
hydrides and hydride cations with relatively small bond energies. 
Thus none of the species S, SH, S^, SH^, or H2S^ can undergo 
an exothermic H atom abstraction reaction with H2. (In other 
words, X-HH2 ^ XH+H is endothermic forX = S, SH, S+, SH+, 
or H2S'^.) Figure 3 illustrates the thermochemistry of the sulphur 
bearing hydrides by displaying the heats of formation of various 
sulphur-bearing species (with associated hydrogen atoms and H2 
molecules involved in their formation.) The endothermicities of 
the reactions X + H2 — > XH + H, expressed in temperature units, 
AE/k, are 9641, 6984, 101 17, 7634, and 4096 K respectively for 
X = S, SH, S+, SH+, and H2S+. 

Clearly, these H atom abstraction reactions leading to the for- 
mation of the sulfur-containing hydrides are negligibly slow at 
the temperatures (< 100 K) typical of diffuse or dense molecular 
clouds. Nevertheless, surprisingly high column densities of H2S 
have previously been observed in both diffuse molecular clouds 
and dense regions of active star-formation. For example, based 
on observations of diff'use foreground material along the sight- 
lines to Sgr B2 (M) and W49N, Tieftrunk et al. (1994) measured 
H2S abundances (of several xlO '' relative to H2 that were con- 
siderably larger than those predicted in standard chemical mod- 
els. They suggested that H2S might have been enhanced by high 
temperature reactions in shocks (e.g. Pineau des Forets, Roueff, 
& Flower 1986) or by hydrogenation reactions on grain surfaces. 
Some smaller H2S abundances (of several times lO"'") were ob- 
served in diffuse clouds at high Galactic latitude by Lucas & 
Liszt (2002), while even larger H2S abundances (up to ~ 10"^) 
have been infeiTed for dense regions of active star-formation 
(Minhetal. 1991). 

More recently, enhanced SH"^ abundances ~ few x 10"^ have 
been measured in the difiiise ISM, using APEX (Menten et al. 
2011) and HerschelfHTPl (Godard et al. 2012). These measure- 
ments have been interpreted as resulting from an increased reac- 
tion rate for S^ -)- H2 — » SH^ -H H in turbulent dissipation regions 
where the gas temperature is elevated and significant ion-neutral 
drift is present. The small SH/H2S abundance ratio ~ 0.13 in- 
ferred in §3 above for the 39kms-' absorbing cloud suggests 
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Fig. 3. Heats of formation of various S-bearing species (with as- 
sociated H atoms produced during their formation), based upon 
thermochemical data from the NIST Chemistry Web book. 



that the reaction SH + H2 -> H2S + H must be similarly in- 
creased, and argues that endothermic neutral-neutral reactions 
are enhanced along with ion-neutral reactions. As described be- 
low, this conclusion is supported by detailed modeling of (1) 
standard photodissociation regions (PDRs); (2) turbulent dissi- 
pation regions (TDRs); and (3) 'C- and 'J' -type shock waves. 

With the use of the Meudon PDR model (Le Petit et al. 
2006), we have found that standard PDR models - i.e. models 
that do not include ion-neutral drift or additional heating mech- 
anisms beyond those associated with ultraviolet irradiation and 
cosmic rays - greatly overpredict the SH/H2S ratio in low den- 
sity diffuse gas, yielding typical values of 10"*, and greatly under- 
predict the SH/H2 and H2S/H2 ratios. We found that SH/H2S 
ratios as small as unity are predicted only in gas clouds with 
densities, nu, larger than lO"* cm~^, and visual extinctions, Ay 
larger than 4, values that are both unreasonably large for the 
39kms"' absorbing cloud toward W49N. We have also inves- 
tigated PDR models that include the effects of grain surface hy- 
drogenation of S and SH and the subsequent photodesorption 
of H2S (Vasyunin & Herbst 2011); these too underpredict the 
SH/H2 and H2S/H2 ratios, unless an unreasonable large gas 
density is assumed. Using models for the chemistry of sulfur- 
bearing molecules in TDRs (Godard et al. 2009) and 'C- and 
'J' -type shocks (Pineau des Forets et al. 1986), we determined 
that the predicted SH/H2S ratios in diffuse molecular clouds are 
~ 10, significantly smaller than those for standard PDRs models 
without grain surface reactions, but nevertheless a factor ~ 100 
larger than what is observed. Here, the production of SH is en- 
hanced by the sequence S-'(H2,H)SH+(H2,H)H2S-'(e,H)SH, in 
which the first two endothermic reactions are significantly en- 
hanced by ion-neutral drift. While the predicted SH and SH"^ 
abundances are broadly consistent with the observations, H2S 
is underpredicted by a factor of 100. As in the PDR models, 
smaller SH/H2S abundances are only predicted in clouds with 
implausibly large tiu and Ay. Apparently, standard models for 
PDRs, TDRs, or 'C- or 'J'-type shocks cannot account for the 



observed SH/H2S ratio, because they fail to predict a sufficient 
enhancement in the rates of neutral-neutral reactions. 

We note however that existing models for sulfur chemistry in 
TDRs and 'C'-type shocks assume that all neutral species share a 
common velocity. As originally pointed out by Flower & Pineau 
des Forets (1998) in the context of CH and CH"^ chemistry, neu- 
tral molecules that are produced by dissociative recombination 
of molecular ions in environments where there is significant ion- 
neutral drift - such as S and SH in the case of present interest 
- initially carry an imprint of the velocity of the ionized parents 
from which they formed. This effect may enhance endothermic 
reactions with H2 that occur before the newly-formed neutral 
species has undergone sufficient elastic collisions to acquire the 
velocity of the neutral fluid. Detailed modeling, which we defer 
to a future pubhcation, will be needed to determine whether the 
SH/H2S ratio serves as a diagnostic of this effect. 
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